Randomized Cellular Automata 

Nino Boccara 
Department of Physics 
University of Illinois at Chicago 
and 

DRECAM-SPEC, CEN Saclay 
91191 Gif-sur-Yvette, France 
boccara@uic . edu 

Abstract 

We define and study a few properties of a class of random automata net- 
works. While regular finite one- dimensional cellular automata are defined 
on periodic lattices, these automata networks, called randomized cellular au- 
tomata, are defined on random directed graphs with constant out-degrees 
and evolve according to cellular automaton rules. For some families of rules, 
a few typical a priori unexpected results are presented. 

1 Introduction 

The purpose of this paper is to define and study a few properties of a class 
of random automata networks. Automata networks are discrete dynamical 
systems which may be defined as follows. Let G = {V, E) be a directed graph, 
where V is the set of vertices and E the set of edges, i.e., a set of ordered 
pairs of vertices. An automata network on \^ is a triple (G, Q, {fi \ i G V}), 
where G is a directed graph on Q a finite set of states usually equal to 
{0, 1, 2, . . . , g — 1}, and fi : Q'^*' t— > Q a mapping called the evolution rule 
associated to vertex i. Ui = {j E V \ {i,j) G E} is the neighborhood of 
i, i.e., the set of vertices connected to i, and \Ui\ is the number of vertices 
belonging to f/j. 
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liV — Zl, the set of integers modulo L, Ui — {i—r^, i—ri+l, . . . , i, i+ 
1, . . . .i + Tj.), i.e., neighborhoods are translation- invariant, and the mappings 
fi are not site-dependent and all equal to /, the corresponding automata 
network is a finite one- dimensional n-input cellular automaton (CA) of size 
L, where and are, respectively, the left and right radii of rule / and 
n — re+rr+1- The global state of such a CA is represented by a configuration, 
i.e., an application c : Q. The space C of all configurations is, therefore, 

equal to Q^^. If s{i,t) denotes the state of vertex i at time t, its state at 
time t -|- 1 is given by 



The apphcation St : i ^ s{i,t) defined on Zl is the state of the CA at time 
t. Since the state St+i is entirely determined by the state St at time t and 
the rule /, there exists a mapping F : C ^— > C such that St+i = F{St). That 
is, if c : i I— > c{i) is a configuration, then its image F{c) is the configuration 



F is referred to as the CA global rule. As time t increases, the set F*(C) 
evolves towards the limit set Ap, i.e.. 



Apis the maximal F-invariant subset of C, i.e., if X C C is such that F{X) = 
X, then X C Ap. 

On the set of all one-dimensional g-state n-input CA rules defined on 
the operators R and C, respectively called reflection and conjugation and 
defined by 

R f{xi, X2, ■ ■ ■ , Xn) = f{Xn, Xn^i, . . . , Xi) 

C f{xi, X2,...,Xn)^f{q-l-Xi,q-l-X2,...,q-l-Xn), 

are the generators of a group of transformations that allows to partition the 
set of rules into equivalence classes. Rules belonging to the same class have 
similar properties. 

The temporal evolution of a one-dimensional CA of size L is represented 
by its spatiotemporal pattern, i.e., the two-dimensional subspace 



s{i,t+ 1) = f(s{i -re,t),s{i -re + l,t), . . . ,s{i + rr,t). 



(1) 



F{c) : i ^ f{c{i - r^), c(i - -M), . . . , c(i r^)). 




{{i,t) |ieZi,ieN}, 
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where N is the set of nonnegative integers. 

The simplest CAs are the so-called elementary cellular automata in which 
the finite set of states is Q = {0, 1} and the rule's radii are both equal to 1. 
There exist 2^ = 256 different elementary CA rules / : {0, 1}'^ — > {0, 1}. The 
rule of an elementary CA can be specified by its look-up table that gives the 
image of each of the eight three- vertex neighborhoods. That is, any sequence 
of eight binary digits specifies an elementary CA rule. Here is an example: 



111 


110 


101 


100 


oil 


010 


001 


000 


1 





1 


1 


1 












Following Wolfram [Ij, a rule code number may be associated with each 
CA rule. If Q = {0, 1}, this code number is the decimal value of the binary 
sequence of images. For instance, the code number of the rule above is 184 
since 

IOIIIOOO2 = 2^ + 2^ + 2^ + 2^ = 184io. 

More generally, the code number N{f) of a one- dimensional | Q|-state n-input 
CA rule / is defined by 

If, for all i G Z^, each vertex of the neighborhood Ui except i is replaced 
with a probability p by a randomly selected vertex in Z^, in such a way 
that the vertices of every resulting neighborhood are different, we call the 
corresponding automata network a randomized cellular automaton (RCA). 
Neighborhoods {Ui \ i E I^l} are no more translation-invariant. An example 
of randomized directed graph is represented in Figure [TJ 

This process of reconnecting vertices to different "neighbors" has been 
used by Watts and Strogatz to build up the first example of small-world 
network [2]. In the case of a RCA there is, however, an important constraint. 
The out-degrees {dout{i) I i ^ Z^} of all the vertices of the randomized 
directed graph on which is defined the new dynamical system are all equal, 
that is, for all i G Z^, dout{i) = \Ui\ — 1 = re + rr- Note that the values of 
din and dout do not take into account the links | « G V}. On the other 

hand, the in-degree of a given vertex is a random variable Di^ that can be 
written as the sum Di + D2 of two independent binomial random variables. 
Di corresponds to the number of original incoming links that have not been 



3 



7 



8 



Figure 1: A randomized directed graph: V = Ziq, n = 3, p = 0.6. The links 
{{i,i) I « G V} have not been represented. 

reconnected to another vertex during the randomization process, and D2 to 
the number of new incoming hnks. The generating function Gin of Di^ is, 
therefore, given by 



As mentioned above, L is the size of the original one-dimensional CA (i.e., 
the number of vertices of the graph) supposed to be large. Hence, the average 
value and the variance of Di„ are 



with 



(p+ (1 -p)z) 




(An) =ri + rr 
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If L is large and p is not too small, i.e., if pL is large, then D2 is well- 
approximated by a Poisson random variable of parameter A = + r^. 

If, after randomization, the neighborhood Ui = (jV^ , • • • ■ ■ ■ , jr^), 

the state s{i, t + 1) of vertex i at time t + 1 is given by 

S{i, t + l) = f{s{jr,,t), . . . , S(j_i, t), S{i, t), S{ji, t), . . . , S{jr,,t)) . (2) 

Translational invariance, which a characteristic property of regular CAs, is 
lost. 

2 Mean-field behavior 

For a given value of the randomizing probability p, statistical averages over 
many different randomized graphs evolving according to a given rule / should 
tend, as p increases, towards values predicted by the mean-field approxima- 
tion applied to the corresponding regular CA. We shall see in the next section 
that this not always the case. In this section we first present a few CA families 
that exhibit this mean-field behavior. 

2.1 Surjective CA rules 

A global CA rule F is surjective if any configuration c G C has a preimage. 
This implies F{C) = C. For a two-state CA, it can be shown [3] that all 
blocks of zeros and ones of a given length have exactly the same number of 
preimages. As a consequence, the average value of the density of ones is equal 
to |. Since this result is correctly predicted by the mean- field approximation, 
the randomization process should have no effect on the value of the average 
density of ones which should be equal to | for all values of p. 

There exist 30 elementary surjective CA rules. Six of them (rules 15, 
51, 85, 170, 204, and 240) are trivial in the sense that their CA rules 
/ : {xi,X2,X3} I— s> f{xi,X2,X3) are functions of one variable only. For the 
remaining 24 rules (rules 30, 45, 60, 75, 86, 89, 90, 101, 102, 105, 106, 120, 
135, 149, 150, 153, 154, 165, 166, 169, 180, 195, 210, and 225) we verify that 
the randomization process, while it destroys the structure of the spatiotem- 
poral pattern, does not change the value of the density of ones, equal to | for 
all values of p. It can also be verified that, for not very large CA size (1000), 
even the statistics of blocks of various lengths also appears to be independent 
of p. 
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2.2 RCA rule 54 



CA rule 54 is defined by the following map 

/54(a;i,a;2,a;3) = 



if (xi, xa, xs) = (0, 0, 1), (0, 1, 0), (1, 0, 0, ) or (1, 0, 1) 
otherwise. 



Rule 54 is one of the most complex elementary CA. Its temporal evolution is 
very slow. It is found that the density of ones tends towards its asymptotic 
value (exactly |) as t~^-^^ [4J. Its spatiotemporal pattern may be viewed 
as particlelike structures (defects) of 4 different types evolving in a regular 
background (Figure [2] ). These particles present a rich variety of interac- 
tions [4j which has received a very nice group-theoretical interpretation due 
to B. Martin [5]. These particles move in a periodic background whose space 
and time periods are both equal to 4. This background can be viewed as a 
square lattice whose elementary cell is the 4x4 tile: 

111 

1 

110 1 ^"^^ 
10 

The definition of the 4x4 tile is not unique. There are actually 8 different 
4x4 tiles according to the 8 possible choices of the first row that are [5]: 

Vo = (0,0,0,1) 6o,i = (1,0,1,1) 

6i,o = (0,0,1,0) = (0,1,1,1) , . 

62,0 = (0,1,0,0) 62,1 = (1,1,1,0) ^ ^ 

63,0 = (1,0,0,0) 63,1 = (1,1,0,1) 

Tile [3] corresponds to background 61^1. As shown by Martin, the different par- 
ticles appear as frontiers between different backgrounds, and can be labeled 
by the elements of the translation group of the regular background. 

If, starting from a configuration belonging to the limit set, we rewire 
one edge (replacing + 1) by with j i + 1), after a few itera- 

tions (whose number depends upon the states of the vertices located in the 
vicinity of vertex j) we may obtain very different results. In particular, as 
shown in Figure [3l we may observe rather complicated structures associ- 
ated to (yf-particles continuously radiating pairs of w-particles (on particles' 
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Figure 2: The four different particles found in the spatiotemporal pattern of 
CA54. 

nomenclature refer to [1]). Note, on the right figure, the existence of a pair 
of new defects appearing after one iteration. 

Increasing the randomizing probabihty p eventually destroys the regular 
spatiotemporal pattern of CA 54. Starting from a random configuration, 
the structure of the spatiotemporal pattern in the infinite time limit clearly 
shows that the asymptotic value of the density of ones is equal |. This 
result is correctly predicted by the mean-field approximation but, since the 
evolution of CA 54 is extremely slow, it takes a very large number of iterations 
to approach this limit. Increasing the probability p, notably increases the 
convergence towards this limit. 

2.3 RCA 18 

For a given value of the randomizing probability p, statistical averages over 
many different randomized graphs evolving according to a given rule / should 
tend, as p increases, towards values predicted by the mean-field approxima- 
tion applied to the corresponding regular CA. This is the case of RCA 18. 
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Figure 3: Pairs of g -particles radiating w -particles. 



Rule 18 is defined by the following map 

fis{xi, X2, X3) = 



if (xi,X2,X3) = (0,0,1) or (1,0,0), 
otherwise. 



The temporal evolution of CA 18 is well understood and rather simple. Con- 
figurations belonging to the limit set consist of sequences of zeros of odd 
lengths separated by isolated ones. The average number of sequences of 
length 2?7. + 1 per vertex is equal to 1/2""^^ [6J. This implies that, in the limit 
of large L, the fraction of ones of a configuration belonging to the limit set 
is exactly equal to ^. In the spatiotemporal pattern a sequence of two ones 
or a sequence of zeros of even length is a defect. Since two sequences of zeros 
of odd lengths separated by two neighboring ones generate, at the next time 
step, a sequence of zeros of even length, configurations generated by rule 
18 contain defects of one type only. During the evolution these defects be- 
have as moving particles. They perform pseudo random walks and annihilate 
pairwise when they meet. This process has been studied by Grassberger [7] 
who found that, starting from a random initial configuration, the density of 
defects decreases as If, following [4|, we view these defects as particle- 

like structures evolving in the background represented by the spatiotemporal 
pattern generated by configurations belonging to the limit set, the pairwise 
annihilation may be viewed as the annihilation of a particle-antiparticle pair, 
the antiparticle, in this case, being identical to the particle. 
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If, starting from a configuration belonging to the limit set, we rewire 
only one edge, that is, we replace, for instance, the edge (i, i + 1) by the edge 
(z, j), where j ^ {i — 1, + after a few iterations (whose number depends 
upon the states of the vertices located in the vicinity of vertex j), a pair of 
defects is created close to site i. As a consequence of the structure of the 
limit set of CA rule 18, it is not possible to create a single defect, only pairs 
of defects may be created. Rewiring randomly more vertices creates more 
defects, some being of a type not found in the spatiotemporal pattern of CA 
18, and eventually the structure of this pattern is completely destroyed. 

The asymptotic value of the fraction of ones varies as a function of the 
randomizing probability p (Figure H] ). For p = 1, which corresponds to 
maximum disorder, the density of ones is close to its mean-field value (0.293). 

0.3 
0.29 

■H 0.28 

to 
c; 

^ 0.27 
0.26 
0.25 

0.2 0.4 0.6 0.8 1 

P 

Figure 4: Densities of ones as a function of the randomizing probability p for 
RCA 18. 



3 Number-conserving rules 

A one-dimensional g-state n-input CA rule / is number- conserving if, for all 
cyclic configurations of length L > ra, it satisfies 

f{x, , X2, . . . , Xn-l, Xn) + f{x2, X3, . . . , Xji, Xji+l ) I ' ' ' 

+ f{xL, Xi . . . , Xn-2, Xn-l) = Xi + X2 H h Xl- (5) 

It can be shown that [8] 
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Theorem A one- dimensional q-state n-input CA rule f is number-con- 
serving if, and only if, for all {xi,X2, ■ ■ ■ , Xn) G Q", it satisfies 



n-l 



f{Xi,X2,..., 




-/( 0,0, .,0^ ,Xi,X2, ■ ■ .,Xn-k)), (6) 

k 

There exist only two number-conserving elementary CA rules, namely 
rule 184 and its conjugate rule 226. The mean-field map of rule 184 is 
p hH> + 2p^(l — p) + p(l — p)^ = p, where p is the density of ones. This 
result, which is valid for all number- conserving rules, could suggest that the 
density of ones should not be affected by the randomization process, i.e., the 
density of ones should, at least in the steady state, remain constant. Our nu- 
merical simulations show that this is not the case: the steady state behavior 
is extremely dependent upon the graph structure. 

For a given value of the randomizing probability p, the average density 
over all vertices becomes either eventually periodic or randomly distributed. 
If p is small, the mean value of the the density remains close to the average 
density of the initial configuration as shown in figure [51 




time time 

Figure 5: RCA 184: average density of ones as a function of time after a 
long transient for a randomizing probability p = 0.05 for two different random 
structures. Graph size: 1000, initial density: 0.5. Eventually periodic density 
with period 6 (left figure ) and randomly distributed density with mean value 
equal to 0.469 and a standard deviation equal to 0.00344 (right figure). 

For larger values of p, when the average density is eventually periodic, 
the period is very sensitive to the graph structure. For example, for p = 1, 
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changing the structure, periods may takes values such as 2, 6, 17, 61, 304 or 
1547, and starting from an initial configuration with an average density of 
ones equal to 0.7, the mean value of the oscillating density can be as low as 
0.13 or as high as 0.97. When the average density is randomly distributed, 
the probability distribution is also extremely sensitive to the graph structure. 
As shown in figure M the probability density is sometimes well approximated 
by a Gaussian distribution, but this is far from being always the case. 



Figure 6: RCA 184: Example of a Gaussian probability distribution of the 
density of ones for p = 1. Lattice size: 1000, density of the initial configura- 
tion: 0.5, mean value: 0.547, standard deviation: 0.0317. 

When the probability p is very small, defects in the spatiotemporal pat- 
tern are distant. In the steady state, their numbers remain constant and 
they interact periodically in time. Hence, the number of active sites is peri- 
odic as a result of the spatial periodicity (Figure [7j). Increasing probability 
p create more defects and their respective distances in a large lattice are 
not commensurate, the resulting average density of active sites is no more 
periodic. 

Similar results are found for 4-input number-conserving rules except for 
rules 51148, 56528, 57580, and 62660 which belong to the same equivalence 
class. These four RCAs exhibit a totally different behavior, they are even- 
tually number- conserving, that is, after a finite number of iterations, for all 
values of the randomizing probability p and all graph structures, the average 
density becomes time-independent. The final average density depends upon 
the value of the average density of the initial configuration, the randomizing 
probability p and, for a given p, upon the graph structure. 
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Figure 7: One- and two-defect RCA184, lattice size: 100. 

4 Conclusion 

We have defined and studied a few properties of a class of random automata 
networks. The spatial structure of these networks is a directed graph with 
vertices having a constant out-degrees and random in-degrees resulting from 
a rewiring process identical to the rewiring process used by Watts and Stro- 
gatz to define the first example of small-world network. When these sys- 
tems evolve according to cellular automaton rules for different values of the 
randomizing probability, we find a rich variety of unexpected behaviors, in 
particular in the case of number-conserving cellular automaton rules. 
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